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ABSTRACT 

We investigate p-process nucleosynthesis in a supercritical accretion disk around a compact object 
of 1.4 Mq, using the self-similar solution of an optically thick advection dominated flow. Supercritical 
accretion is expected to occur in a supernova with fallback material accreting onto a new-born compact 
object. It is found that appreciable amounts of p-nuclei are synthesized via the p-process in supernova- 
driven supercritical accretion disks (SSADs) when the accretion rate m = Mc^/(16LEdd) > 10®, where 
^Edd is the Eddington luminosity. Abundance prohles of p-nuclei ejected from SSADs have similar 
feature to those of the oxygen/neon layers in Type II supernovae when the abundance of the fallback gas 
far from the compact object is that of the oxygen/neon layers in the progenitor. The overall abundance 
prohle is in agreement with that of the solar system. Some p-nuclei, such as Mo, Ru, Sn, and La, are 
underproduced in the SSADs as in Type II supernovae. If the fallback gas is mixed with a small fraction 
of proton through Rayleigh-Taylor instability during the explosion, significant amounts of ®^Mo are 
produced inside the SSADs. ®®Ru and ^®®La are also produced when the fallback gas contains abundant 
proton though the overall abundance profile of p-nuclei is rather different from that of the solar system. 
The p-process nucleosynthesis in SSADs contributes to chemical evolution of p-nuclei, in particular ®^Mo, 
if several percents of fallback matter are ejected via jets and/or winds. 

Subject headings: Accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances — stars: 
supernovae: general 


1. INTRODUCTION 

There exist 35 neutron dehcient stable nuclei with mass 
number A > 74 referred to as p-nuclei. These nuclei 
are found in meteorites only inside the solar system; iso¬ 
topic anomalies involving these nuclei are also found in 
some primitive meteorites. Many production sites of 
the p-nuclei have been proposed: oxygen/neon layers of 
highly evolved massive stars during their presupernova 
phase (Arnould 1976) and during their supernova explo¬ 
sion (Woosley & Howard 1978); X-ray novae (Schatz et 
al. 1998); neutrino-driven winds originated from a nascent 
neutron star shortly after supernova explosion (Hoffman 
et al. 1996); Type la supernova explosion (Howard, Meyer 
& Woosley 1991); helium-accreting CO white dwarfs of 
sub-Chandrasekhar mass (Goriely et al. 2002). Among 
them, the most promising is the oxygen/neon layers during 
Type H supernova explosion. The p-nuclei are synthesized 
by the photodisintegrations of s-nuclei (s-process seeds) 
produced in the layers during the core helium burning in 
the progenitor. Photodisintegrate (7, n) reactions are fol¬ 
lowed by (7,p) and/or (7,0;) reactions. The production 
of p-nuclei via the subsequent photodisintegrations is re¬ 
ferred to as a p-process. Extensive investigations of the p- 
process in Type H supernovae are performed by Prantzos 
et al. (1990b), Rayet, Prantzos & Arnould (1990), Rayet et 


al. (1995), and Arnould, Rayet & Hashimoto (1998) with 
a large nuclear network. A p-process model of Rayet et 
al. (1995) in the oxygen/neon layers during Type H super¬ 
novae reproduces the overall abundance profile of p-nuclei 
in the solar system. However, the model has a conspicu¬ 
ous shortcoming of underproduction of ®^’®^Mo, ®®’®®Ru, 
and ^®®La. Abundant productions of such p-nuclei have 
been invoked to neutrino processes (Woosley et al. 1990; 
Goriely et al. 2001); Yield of ^®®La in the oxygen/neon 
layers is increased via neutrino interactions though those 
of Mo and Ru are not appreciably changed (Goriely et 
al. 2001). Enhancement of s-process seeds synthesized in 
the layers during core helium burning could improve the 
underproduction of Mo and Ru with an extremely high 
rate of ^^Na(a, n)^®Mg reaction (Costa et al. 2000). How¬ 
ever, such a high rate is possibly excluded by a recent 
experiment (Jaeger et al. 2001). 

During supernova explosion in massive stars, certain 
fraction of ejected matter is expected to fall back onto a 
new-born compact object with supercritical accretion rates 
greater than lO^Mcr (Woosley & Weaver 1995, hereafter 
WW95), where Mci is the critical mass accretion rate given 
by 16LEdd/c^ with the Eddington luminosity, LEdd- If 
the fallback material has substantial angular momentum, 
formation of an accretion disk is inevitable (Mineshige et 
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al. 1997). Fujimoto et al. (2001) have investigated nu¬ 
cleosynthesis inside supernova-driven supercritical accre¬ 
tion disks (SSADs) and shown that appreciable nuclear 
reactions take place if the accretion rate M > lO^Mcr- 
At such high accretion rates, SSADs contain the regions 
where density, temperature, and dynamical time-scale are 
comparable to those of the oxygen/neon layers in Type II 
supernovae. Therefore, the p-process is likely to operate 
efficiently inside SSADs. 

We propose a SSAD as a candidate for the produc¬ 
tion of p-nuclei. In §2, we present a model of SSADs 
and a nuclear network to calculate the abundance evo¬ 
lution inside a SSAD. We also discuss possible progenitors 
for p-process nucleosynthesis inside SSADs and the initial 
chemical abundances of the fallback material far from the 
central compact remnant. In §3, we describe abundance 
distributions of p-nuclei inside SSADs, ejected masses of 
these nuclei from SSADs, and effects of proton and helium 
contamination of the fallback gas via Rayleigh-Taylor in¬ 
stability. We discuss our results in §4 and finally present 
concluding remarks in §5. 

2. MODEL AND INPUT PHYSICS 

2.1. The Supercritical Accretion Disk Model 

Physical quantities of supercritical accretion disks with 
M ^ Mcr are well described in terms of the self-similar so¬ 
lution. The density p, temperature T and drift time-scale 
tdr of accreting gas are evaluated with the self-similar disk 
model as adopted by Fujimoto et al. (2001) to be 



where r is the radial coordinate, M is the mass of the 
central compact object, Ovis is the viscous parameter, 
TO = M/Mcn and rg = 2GM/c^ is the Schwarzschild ra¬ 
dius. The model has been assumed to be steady and ax- 
isymmetric with the self-gravity ignored. It is noted that 
in SSADs with toSIO^^ advective cooling dominates over 
radiative and neutrino cooling; it balances with viscous 
energy generation. A disk is also assumed to be extended 
from the last stable circular orbit around a black hole to 
the location of accretion shock ~ 10®cm (Fryer, Colgate & 
Pinto 1999; MacFadyen, Woosley & Heger 2001). Hence, 
the inner and outer boundaries of the disk, rin and rout) 
are set to be 3 rg and 3x 10^ rg (= 1.2 xIO^M/Mq cm), re¬ 
spectively. A sequence of the models is specified in terms 


of M, Ovis) and to. In the present paper, we fix both 
M = 1.4 Mq and Ovis = 0.01 (e.g, Hawley 2000), while to 
is left to be a parameter. 

2.2. The Nuclear Reaction Network 

We have developed a nuclear reaction network which 
has been extended from a smaller network (Hashimoto 
& Aral 1985; Koike et al. 1999). The nuclear data are 
taken from the data base REACLIB^, which includes many 
experimental rates and theoretical Hauser-Feshbach rates 
(Rauscher & Thielemann 2000, 2001). We have also used 
the experimental mass data (Audi & Wapstra 1995). The 
reactions included in the network are ( 0 , 7 ), (n, p), ( 0 , 0 ), 
(P)7)) (P)<t), ( 0 , 7 ), 30, and their inverse reactions. The 
various channels of ^®C-|-^®C, ^^C-|-^® 0 , and ^® 0-|-^®0 are 
also included in the network. For and decays, REA- 
CLIB rates are updated with the corresponding data in 
Horiguchi et al. (1996) if available. We have also added 
and /S®" decays and electron and positron captures (Fuller 
et al. 1980, 1982a, 1982b). Some rates are revised with 
experimental rates (see details in Koike et al. 1999). Our 
network includes 1988 nuclides from neutron and proton 
up to ®®®Bi, presented in Table 1. The network is solved 
implicitly with an inverter of a sparse matrix (Timmes 
1999). 

2.3. Possible Progenitors for P-Process Nucleosynthesis 
inside SSADs 

As the case of Type H supernovae, synthesis of p-nuclei 
inside SSADs requires the s-process seeds. The fallback 
after a supernova explosion is likely to be induced by deep 
gravitational potential of relatively massive stars and/or 
reverse shock inwardly propagating from the outer compo¬ 
sition interfaces (Colgate 1971; Chevalier 1989; WW95). 
A total amount of fallback material increases as mass of 
a progenitor are enhanced (WW95; Fryer 1999, hereafter 
F99; Fryer & Heger 2000); Massive fallback is likely to take 
place for the progenitors with Mms ~20Mq due to their 
large binding energy of the iron core, where Mms is the 
mass of the progenitor on the main-sequence. For massive 
stars with Mms~ 4 OM 0 , supernova explosion cannot take 
place and collapsed matter forms a black hole promptly 
(F99). Several seconds after the core collapse a quasi¬ 
steady accretion disk is formed around the hole if the pro¬ 
genitor has sufficient angular momentum (MacFadyen & 
Woosley 1999), which is in a reasonable range of a progeni¬ 
tor model of Heger, Langer & Woosley (2000). Highly vari¬ 
able supercritical accretion with 0 . 01 - 0.1 Mq s“^ is main¬ 
tained for approximately 10-20 s (MacFadyen & Woosley 
1999). For such high accretion rates, however, the accre¬ 
tion disk has enough high temperature (710^® K) for ac¬ 
creting gas to contain only protons and neutrons (Popham, 
Fryer & Woosley 1999). Therefore, the progenitors with 
Mms ~ 4 OM 0 is irrelevant for p-process nucleosynthesis in¬ 
side SSADs. 

While for progenitors with Mms = 2 O- 4 OM 0 a neu¬ 
tron star is formed soon after the core collapse. Super¬ 
nova launches outward-going shock successfully but some 
fraction of supernova ejecta falls back onto the nascent 
neutron star, which is transformed to a black hole via 
massive fallback (F99). The accreting matter onto the 


The REACLIB compilation of reaction rates is available at: ftp://quasar.physik.unibas.ch/pub/tommy/astro/reaclib/ 
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hole forms a disk with 1 O“‘*- 1 O“^M 0 s“^ lasting for hun¬ 
dreds to thousands of seconds (MacFadyen et al. 2001). 
After most of silicon layers of the progenitors has fallen 
back into the hole, the accretion rate is declined as 
(Chevalier 1989; MacFadyen et al. 2001) and compositions 
of the accreting matter are those of explosively burned 
oxygen/neon layers via supernova shock. In a 20Mq pro¬ 
genitor model (Hashimoto 1995), the position of rout (= 
SOOOcg) is far from the core and located at Mr — 2 .IM 0 , 
or the oxygen/neon layer, where the peak temperature 
during the shock propagation is too low to proceed ap¬ 
preciable nuclear burning of s-process seeds (S2 x 10® K). 
The material still leaves the s-process seeds even after the 
supernova explosion. Appreciable amounts of the seeds 
can fall back onto the compact object through the accre¬ 
tion disk. Accordingly, p-process nucleosynthesis operates 
inside a supercritical accretion disk driven by supernova 
explosion of the progenitors with Mms = 2 O- 4 OM 0 . 

2.4. Initial Abundances of Fallback Matter 

Chemical composition of fallback matter depends on 
masses of progenitors and fallback mechanisms. Because 
p-process nucleosynthesis inside SSADs needs s-process 
seeds, we set the compositions of the accreting gas at rout 
to be averaged compositions of oxygen/neon layers of the 
progenitors of 10, 20, 30, and 40 Mq, whose models are 
refererred to as MIO, M20, M30, and M40, respectively. 
These abundances were calculated with a network which 
contains 440 nuclei up to Bi. (Prantzos, Hashimoto & 
Nomoto 1990a). Although the progenitors with masses of 
Mms < 2 OM 0 are unlikely to be relevant for p-process 
inside a SSAD, we also calculate abundances of p-nuclei 
for MIO because the abundance pattern is different from 
the others (Figure 1). It should be emphasized that these 
abundances were also adopted as the seed abundances in 
the p-process models of Rayet et al. (1990, 1995). Figure 
1 shows the initial abundances of MIO, M20, and M40, 
which are denoted by the dashed, solid, and thick-solid 
lines, respectively. The initial abundances of M30 are in 
between those of M20 and M40. The solar abundances are 
also depicted with the filled circles (Anders & Grevesse 
1989). 

Compositions of fallback matter are not those of the 
oxygen/neon layers of the progenitors but those of explo¬ 
sively burned layers because of supernova shock heating. 
Abundances of p-nuclei produced inside SSADs depend on 
amounts of s-process seeds, as we can see later. The ex¬ 
plosively burned layers where the peak temperature of a 
shock wave is lower than 2xlO®K still contain compara¬ 
ble s-process seeds to those before the explosion. Accord¬ 
ingly, the initial compositions of the oxygen/neon layers of 
the progenitors are relevant for our p-process model inside 
SSADs at least for fallback of the layers with relatively low 
peak temperature < 2 x 10®K. 

The compositions of fallback material discussed above 
are possibly altered when large-scale matter mixing is sig¬ 
nificant in supernova explosion. Such mixing is probably 
associated with fallback. Hence, we shall discuss effects of 
matter mixing on p-process nucleosynthesis in §3.3. 

3. NUCLEOSYNTHESIS INSIDE SSADS 
3.1. Abundance Distributions of P-Nuclei inside SSADs 


As fallback material accretes onto a central object, tem¬ 
perature and density of the gas increase (see equations 
( 1 ) and ( 2 )) and thus abundances of the accreting gas are 
changed through a sequence of nuclear reactions from an 
initial composition described in §2.4. Using the nuclear re¬ 
action network and the self-similar solution, we can follow 
the evolution of the chemical composition inside SSADs 
during the infall toward the compact object through post¬ 
processing calculations for MIO, M20, M30, and M40. The 
calculations are carried out from rout to rin. It should be 
noted that the nuclear energy generation is assumed to be 
much smaller than the viscous heating. 

For disks with mSlO®, the maximum temperature of 
disks is less than 2 x 10® K which is the minimum temper¬ 
ature for the p-process to proceed in the so called p-process 
layers (PPLs) in Type H supernovae (Rayet et al. 1990, 
1995). Therefore, p-process nucleosynthesis does not oper¬ 
ate in the disks for such low accretion rates. For higher ac¬ 
cretion rates, the accreting gas has higher temperature and 
the p-process can proceed significantly. The abundance 
distributions of p-nuclei are similar for the disks with var¬ 
ious accretion rates, while the production site of p-nuclei 
shifts toward the outer part of the disk as the accretion 
rates increase. Thus, we show the results of m = 10® as a 
representative case. 

Figure 2 shows the abundance profiles, normalized with 
the solar abundances, of representative p-nuclei for M20. 
The abscissa is the radius of the fallback disk in units of 
Tg. The solid, thick-solid, dashed, thick-dashed, dotted, 
and thick-dotted lines indicate the abundances of ^"‘Se, 
®^Mo, ^®®La, ^®®Gd, ^®®Ta, and ^®®Hg, respectively. At 
r > 100 Tg, temperatures of the accreting gas are too low 
to proceed appreciable synthesis of p-nuclei. The abun¬ 
dances of ^®^Gd and ^®®Ta are the same as the initial val¬ 
ues. At r ~ 100 Tg (T ~ 1.3 X 10® K), ^®®Ta increases 
via neutron capture of ^^®Ta, contained in the fallback 
material initially. As the gas moves to the inner part of 
the disk, the temperature increases and the heavy p-nuclei 
{A > 150), such as ^®^Gd, ^®®Ta, and i®®Hg are rapidly de¬ 
stroyed through ( 7 ,n) reactions. At r^50 rg (T>2.0 x 10® 
K), however, these isotopes are drastically produced. This 
is due to the photodisintegrations of the neutron-rich seed 
elements when the temperature of the accreting gas at¬ 
tains to 2 — 3 X 10® K comparable to PPLs (Rayet et 
al. 1995). Thus, the p-processes operate efficiently inside 
the accretion disk as in the PPLs. The intermediate p- 
nuclei (100 < A < 150) also enhance via the p-process. It 
should be noted that tdr — 0 . 1-1 s in this region is compa¬ 
rable to the time-scale of shock propagation through PPLs 
(Rayet et al. 1995). 

For a region r<30rg (T>2.8 x 10® K), ( 7 ,p) and ( 7 , 0 ;) 
reactions of the heavy and intermediate isotopes dominate 
over ( 7 , n) reactions; while for the light p-nuclei {A < 100), 
such as ^'^Se and ®^Mo, ( 7 , n) reactions are still dominant 
processes. The light p-nuclei are eventually depleted via 
( 7 , p) and ( 7 , a) reactions near 20 Cg (T ~ 3.6 x 10® K). It 
is noted that the abundance profiles for all models are simi¬ 
lar, while the more p-nuclei are produced the more massive 
are progenitors because of the richness of s-process seeds 
(Figure 1). 

3.2. Ejected Masses of P-Nuclei from SSADs 
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Many astrophysical objects involving an accretion disk 
are observed to produce jets and/or winds (e.g., Livio 1999 
and references therein). Hence, as the accreting gas falls 
back onto the central object, some fraction of gas is pos¬ 
sibly ejected from the disk via jets and/or winds. The 
processed material ranging from ri„ to Tej is supposed to 
be ejected. Then we can estimate the averaged mass frac¬ 
tion Xi of the i-th p-nucleus ejected from the disks, 

\ X,(r)Y.(r)rdr, (4) 

where S(r) is the column density of the disk. Here M{rej) 
is the mass of the disk from rin to Tej and the fraction 
of the ejected gas to the accreting gas is assumed to be 
constant in radius. To compare the p-nuclei abundances 
ejected from the fallback disks with the solar abundances, 
we calculate the overproduction factors (OPFs) Fi of 35 
p-nuclei (Rayet et al. 1995): 

F, = (X,/X,,0)/Fo, (5) 

where Xi^ q is the mass fraction in the solar system and 
Fq = 0)/35. In many accretion powered sys¬ 

tems, matter ejection via jets and/or winds could be origi¬ 
nated from the inner region of the accretion disks. There¬ 
fore, we set Tej to be lOOrg (Junor et al. 1999). The chem¬ 
ical composition of the ejected material is assumed to be 
constant throughout a jet and/or a wind. Abundances of 
p-nuclei increase due to /3-decays of their unstable parent 
nuclei after the freeze out. However, abundances of most 
p-nuclei are not appreciably enhanced after the /3-decays 
(see §4.2). 

Figures 3 and 4 show the OPFs of 35 p-nuclei for M20 
and MIO, respectively. It is emphasized that the abun¬ 
dance profiles have similar features to those in Type H su¬ 
pernovae (Rayet et al. 1995). The light p-nuclei, 
and ®®’®®Ru are underproduced in our calculations. This 
is due to the deficiency of their s-process seeds (Fuller 
& Meyer 1995). The intermediate mass p-nuclei, ^^^In, 
^^®Sn, and ^^®La, cannot be produced appreciably for MIO 
and M20. It is noted that ^^^In and ^^^Sn are largely 
enhanced through /3-decays after ejection via jets and/or 
winds (see §4.2). ^^^Gd is underproduced in our calcula¬ 
tions, while a large amount of ^®^Gd is synthesized in the 
stellar s-process (Prantzos et al. 1990a). ^®^Er is also un¬ 
derproduced in our calculations. The overproductions of 
^■^Se (M20) and ^®'^Os (MIO) are attributed to the ample¬ 
ness of their s-process seeds. The OPF of ^^°Te is largely 
increased for MIO compared with that of M20. The OPFs 
of Se, Kr, and Sr for MIO are smaller than those for M20 
but those of p-nuclei with A > 168 are larger. Such pro¬ 
files of the OPFs are also responsible for initial abundance 
distributions of s-process seeds. The OPFs for M30 and 
M40 are similar profiles to those for M20. However, Fq for 
M30 and M40 are larger than Fq for M20 (Table 2). This 
is because more massive progenitor has larger amounts of 
s-process seeds (Figure 1). 

Profiles of the OPFs ejected from SSADs depend on 
initial abundance distributions of s-process seeds. The s- 
process seeds are mainly synthesized during a helium core 
burning stage in massive stars. The s-process nucleosyn¬ 
thesis in massive stars has been extensively investigated 
by many authors (e.g., Prantzos et al. 1990a; Kappeler et 
al. 1994; Rayet & Hashimoto 2000); Uncertainties in some 


key nuclear reactions still preclude precise estimation of 
s-process yields. 

3.3. Proton and Helium Contamination through Matter 

Mixing 

Several observations in SN 1987A support the occur¬ 
rence of large-scale mixing in the ejecta. Newly synthe¬ 
sized ®®Ni near the collapsed core is likely to be mixed up 
into the hydrogen envelope as indicated from observations 
in early phase; X-ray and y-ray light curves (e.g., Kumagai 
et al. 1989 and references therein) and optical spectroscopy 
(e.g., Mitchell et al. 2001 and references therein). Two- 
dimensional simulations of supernova explosions (Muller, 
Fryxell & Arnett 1991; Herant & Benz 1992; Nagataki, 
Shimizu & Sato 1998) indicate that the large-scale mixing 
is caused via Rayleigh-Taylor instability. A recent high- 
resolution simulation by Kifonidis et al. (2000) shows that 
a large amount of ^He is mixed down near the iron core by 
the Rayleigh-Taylor instability. Hydrogen is also shown to 
be mixed down to the core (Hachisu et al. 1990). More¬ 
over, several analyses of observations in SN 1987A are also 
preferable to mixing of the outer hydrogen-rich envelope 
down to the inner metal-rich core; the width of plateau-like 
peak of the bolometric light curves (Shigeyama & Nomoto 
1990), line profiles in late spectral evolution (Kozma & 
Fransson 1998), and the light curve for the hrst 4 months 
(Blinnikov et al. 2000). Thus, appreciable amounts of pro¬ 
ton and helium are likely to be mixed down to the inner 
core via the Rayleigh-Taylor instability in SN 1987A. Ger- 
tain fractions of proton and helium are accordingly mixed 
into the oxygen/neon layers. 

Such mixing is probably generic in core-collapsed su¬ 
pernovae; recent X-ray observations in Gas A (Douvion et 
al. 1999; Hughes et al. 2000), emission line profiles of SN 
1988A (Spyromilio 1991) and SN 1993J (Spyromilio 1994), 
and optical and infrared spectroscopies of SN 1995V (Fas- 
sia et al. 1998) and SN 1998S (Fassia et al. 2001). 

Moreover, recent investigation for anomalous composi¬ 
tion in a companion of the black hole binary Nova Sco sug¬ 
gests that the companion is polluted by material ejected 
in the supernova accompanied with the formation of the 
hole (Israelian et al. 1999). The progenitor in this system 
is likely to experience substantial fallback (~ a few Mq) 
and matter mixing in supernova ejecta (Podsiadlowski et 
al. 2002). Such fallback with mixing is also preferable to 
the large abundance of Zn observed in the very metal-poor 
stars (Umeda & Nomoto 2002). Therefore, the fallback 
gas is expected to be contaminated with some fractions of 
proton and helium via the large-scale mixing. 

To investigate effects of proton and helium contamina¬ 
tion of the fallback material on the p-process nucleosyn¬ 
thesis, we add an arbitrary fraction of proton or helium to 
the initial abundance of M20 and calculate the abundance 
distributions inside SSADs. The initial proton abundance 
of the fallback matter highly depends on large-scale mixing 
and fallback mechanism. There exists no reliable estima¬ 
tion of the initial proton abundance. Analyses of observa¬ 
tions in SN 1987A indicate that the mass fraction of proton 
is 0.01 (Blinnikov et al. 2000) or up to 0.3 (Shigeyama & 
Nomoto 1990). Hence we examine the cases of the initial 
mass fraction of proton Xp = 0.01, 0.1, and 0.3, which are 
referred as models M20P001, M20P01, and M20P03, re- 
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spectively. In these models, p-nuclei are produced via not 
only photodisintegration of s-process seeds but radiative 
proton capture of lighter nuclei. 

Figure 5 shows the OPFs for M20P001 with m = 10®. 
The most remarkable is an enhancement of ®^Mo. Even if 
the fallback material contains a small fraction of proton, 
radiative proton captures as well as photodisintegrations 
are efficient. The light p-nuclei, ^"^Se, ^®Kr, and ®‘*Sr, are 
increased, while the light p-nuclei underproduced in M20, 
such as ®^Mo, ®®Ru, and ®®Ru, are still deficient. The 
intermediate mass p-nuclei, ^^"‘Sn, ^^®Sn, and ^®®La are 
underproduced as M20. The OPFs of the heavy p-nuclei 
are decreased compared with M20 though abundances of 
these p-nuclei are not appreciably changed. This is be¬ 
cause the averaged OPF, Fq, is increased (equation (5)). 
We also find that profiles of the OPFs weakly depend on 
Xp for TO = 10®: the OPFs for M20P01 and M20P03 are 
not significantly changed from those for M20P001. The 
averaged OPFs, Fq, of M20P001, M20P01, and M20P03 
are larger than that of M20 and listed in Table 2. For 
models with initial abundances of MIO, M30, and M40 
added by proton of Xp = 0.01, 0.1, and 0.3, the OPFs 
are similar profiles to model with the initial abundance of 
20Mq model added by the same fraction of proton (i.e., 
same Xp) though the averaged OPFs, Fq, are significantly 
increased for models with massive progenitor (Table 2). 

For TO = 10®, the profiles of the OPFs change as pro¬ 
tons are injected. Appreciable amounts of the p-nuclei 
which are deficient in model M20 are produced in SSADs. 
For Xp = 0.01 (M20P001), ®^Mo and ^®®La are overpro¬ 
duced. The overproductions of ®^Mo and ^®®La are more 
prominent in M20P01 (Figure 6) compared with M20P001. 
When Xp = 0.3 (model M20P03; Figure 7), 9®Ru, ^^^Sn, 
and ^^^Sn become abundant because of very efficient pro¬ 
ton captures. However, the overall profile of the OPFs be¬ 
comes worse compared with the solar one. Fq drastically 
increases compared with the cases of small proton inclu¬ 
sion; Fq are equal to 67.68, 41.84, and 2882, for M20P001, 
M20P01, and M20P03, respectively. It should be noted 
that ®‘*Mo, ®®Ru, and ^^®Sn are underproduced even if 
protons exist abundantly. 

Concerning helium contamination of fallback material 
via large-scale mixing, the OPFs are not appreciably 
changed from those of M20 even if a significant fraction 
of ^He, up to 0.5, is added to the initial abundance of 
M20. Hence, amounts of p-nuclei produced inside SSADs 
are independent of an amount of ^He included in fallback 
material, while a elements, such as ^^Ti, are enhanced. 


4. DISCUSSION 

4.1. Contribution of SSADs to Chemical Evolution of 

P-Nuclei 

We discuss contribution from SSADs to chemical evo¬ 
lution of p-nuclei to estimate relative contribution from 
SSADs in respect to supernovae. In this aim we average 
the overproduction factor Fq of p-nuclei produced by su¬ 
pernovae or SSADs using the initial mass function (IMF) 
f{M). The IMF-weighted OPF is calculated as 


(Fb)lMF 


Eo(M)Meont(M)g(M)dM 

/^“^(M)dM 


( 6 ) 


where M\ and Mu are the lower and upper mass limits of 
progenitors and Mcont is the total mass (in units of Mq) 
which contributes to the production of p-nuclei via super¬ 
novae or SSADs. 

We assume that the progenitors of 13 < Mms/Mq < 20 
can produce p-nuclei via supernova explosion (Rayet et 
al. 1995), while the progenitor with 20 < Mms/Mq < 40 
can yield p-nuclei through SSADs driven by massive fall¬ 
back. For supernovae and SSADs, Mcont is the mass of the 
PPLs, MppL (Rayet et al. 1995) and the ejected mass from 
a disk via jets and/or winds, Mej, respectively. Adopting 
f{M) oc M“^® proposed by Salpeter (1955), we can cal¬ 
culate (Fo)imf as 18.3 for supernovae. Here we use the 
results of Rayet et al. (1995) for Mppp and Fq{M) (in 
their Tables 2 and 3). A calculation is performed in a 
similar manner to Rayet et al. (1995). 

Next we estimate (Fo)imf for SSADs for the two cases: 
no proton mixing and complete mixing. In actual situa¬ 
tions, uncertain fractions of fallback matter are affected 
by matter mixing. Hence (Fo)imf for SSADs is located 
between the two estimations. Firstly we discuss the case 
in absence of matter mixing and then those with complete 
mixing. 

4.1.1. The Case in The Absence of Matter Mixing 

To calculate (Fo)imf for SSADs, we need an estimation 
of the total fallback mass with s-process seeds, Mgeed- An 
amount of fallback material has been estimated by WW95 
and F99 through numerical simulations of core collapse 
and supernova explosion. Amounts of fallback material 
Mf are presented in Table 3 for the fallback models of 
WW95 and F99. To operate the p-process in SSADs, fall¬ 
back material should contain the s-process seeds. Accord¬ 
ingly, to calculate Mgeed, besides the amount of fallback 
matter, both the abundance distribution of the progeni¬ 
tor and the peak temperature of the supernova shock are 
needed. In Table 3, we also show the mass coordinate of 
the base of oxygen burning layers, Mob and the mass co¬ 
ordinate where the peak temperature is 2 x 10® K, Mdest 
for the progenitor and explosion models with masses of 15, 
20, 25, and 40 Mq (Hashimoto 1995). All s-process seeds 
are assumed to be destroyed in the region of the mass co¬ 
ordinate Mr < Mdest and thus disk p-process can operate 
only for fallback material with Mr > Mdest- Hence we 
can estimate Mgeed, which are also presented in Table 3, 
for WW95 and F99. For the fallback model of WW95, 
Mseed(3OM0) is estimated as 2.41M0 provided that Mdest 
for the 30Mq progenitor is equal to that of the 25Mq 
progenitor. For the fallback model of F99, we set sim¬ 
ply Mseed(2OM0) = O.5[Mseed(15M0)-|-Mgeed(25M0)] and 
Mgeed(3OM0) = i[2Mgeed(25M0) +Mgeed(4OM0)]. 

Moreover, some fractions of synthesized p-nuclei inside 
SSADs are ejected from the disks via jets and/or winds as 
in many objects associated with an accretion disk (e.g., 
Livio 1999). We assume that the ratio e of the mass 
ejection rate Mej to the mass accretion rate is constant 
in time, or Mej = (M. Hence the total ejected mass is 
also expressed as Mej = eMseed- The ratio e are 0.001- 
0.01 (Eggum, Coroniti & Katz 1988) and 0.03-0.1 (Ku- 
doh, Matsumoto & Shibata 1998) estimated from numer¬ 
ical simulations of jets originated from an accretion disk. 
However, e may be larger as suggested by observations in 



6 


Fujimoto et al. 


X-ray binaries; SS433 (Kotani 1997) and GRS 1915-1-105 
(Fender & Pooley 2000; Belloni et al. 2000). 

Using the calculated values of Fq for our SSAD models 
(Table 2), we evaluate (Fo)imf as 1.640 eo.oi and 2.864 
Co .01 with the estimation of amounts of fallback matter of 
WW95 and F99, respectively, where eo.oi = e/0.01. Thus 
(Fb)iMF for SSADs are smaller than that for supernovae 
(= 18.30) unless e is as large as 0.1. It should be noted 
that (Fo)imf slightly decreases for a steep IMF slope with 
^(M) (X (e.g., Kroupta 2002) and also reduce by 

30-50% if Mu = 3OM0. 

4.1.2. The Case with Complete Matter Mixing 

We proceed the case of fallback with complete proton 
mixing. We assume matter contains proton of Ap = 0.01 
throughout fallback. The proton mixing needs the hydro¬ 
gen envelope of progenitors, which is decreased due to mass 
loss for massive progenitors (e.g., Chiosi & Maeder 1986). 
The upper limit Mu possibly decreases for cases with pro¬ 
ton mixing and are taken to be Mu = 35Mq (Fryer & 
Kalogera 2001). As in the absence of matter mixing, we 
calculate (Fo)imf as 5.027 eo.oi and 9.583 eo.oi for the fall¬ 
back models of WW95 and F99, respectively. Here we 
set Fo(35Mo) = 0.5[Fo(30Mo) -t Fo(4OM0)]. The en¬ 
hancement of (Fb)iMF compared with no proton mixing is 
responsible for the increase in Fq for the models with pro¬ 
ton mixing (Table 2). Hence, in the case with complete 
mixing, the SSADs-origin p-nuclei are likely to compara¬ 
bly contribute to the chemical evolution of p-nuclei as the 
supernovae-origin p-nuclei. 

In particular, SSADs probably have significant contri¬ 
bution to the ®^Mo evolution. Overproduction of ®^Mo 
takes place for any m (Figures 5-7). Accordingly, the over¬ 
production is realized inside SSADs throughout the dura¬ 
tion in which a large fraction of fallback matter accretes 
onto a central compact object. Even for the cases of small 
Ap, ®^Mo is overproduced inside SSADs; This is the case 
of Ap >0.005. Thus, contribution to the ®^Mo evolution 
is considerable for a small amount of proton ~ O.OIM© 
mixed into the fallback matter of 0.1 to several Mq. 

On the other hand, ®®Ru and ^^®La are overproduced 
only for the case of m = 10®. The accretion rate of fall¬ 
back matter is declined to relatively small accretion rates 
m ~ 10® at several months after the explosion (Mineshige 
et al. 1997) and large amounts of fallback matter have al¬ 
ready accreted onto a central remnant. Therefore, only 
a small amount of fallback matter are processed inside 
SSADs, which are enable to overproduce ®®Ru and ^®®La. 

4.2. Abundance Change of P-Nuclei through Jets and 

Winds 

In SSADs, p-nuclei are mainly produced in the region 
20-70 Tg for m = 10® (Figure 2) where an accreting gas is 
mainly ejected via winds. Parent nuclei which decay into 
p-nuclei through /3-decays are also expelled via winds and 
thus enhance amounts of p-nuclei. Taking into account 
the /3-decays after freeze out, we estimate abundances of 
p-nuclei ejected from a SSAD. We find that p-nuclei other 
than ®^Mo, ®^Mo, ^^®In, and ^^®Sn are not appreciably 
increased via /3-decays of their parent nuclei; Their ratio 
after to before /3-decays are 1.100, 1.215, 13.77, and 18.33, 
respectively, for M20. Here all the parent nuclei has been 


assumed to decay to their daughter p-nuclei. When pro¬ 
ton is abundantly mixed into fallback material, the abun¬ 
dances of ®^Mo and ®^Mo are not significantly changed 
but ^^®In and ^^®Sn are largely enhanced. For M20P03 
(Ap = 0.3), the ratio of the abundances of ^^®In, and 
^^®Sn after to before /3-decays become 9.797 and 42.41, 
respectively. However, the averaged OPFs, Fq, are not 
appreciably changed because of small OPFs of ^^®In and 
^^®Sn up to ~ 0.07 before /3-decays (Figures 3-7). Thus, /3- 
decays after freeze out has only minor effect on estimation 
of (Fb)iMF discussed in §4.1. 

4.3. Neutron Star Kick 

If a compact object breaks out a supernova remnant, 
matter fallback onto the object ceases. Many observations 
reveal that pulsars have larger space velocities (200-1000 
kms“^) than those of their progenitor (e.g., Lai et al. 2001 
and references therein). Such high velocities are invoked to 
kicks of nascent neutron stars via asymmetric supernova 
explosions. The supercritical accretion (to > 10®) investi¬ 
gated in the present paper is realized until about 4 years 
after an explosion (Mineshige et al. 1997). The expansion 
velocity of the supernova remnant is much larger than the 
typical kick velocity of pulsars. The material therefore 
continues to fall back during the supercritical accretion of 
our interest. However, the later expansion of the remnant 
is decelerated by ambient material from progenitor winds. 
The pulsar eventually breaks through the supernova rem¬ 
nant and matter fallback onto it ceases. 

5. CONCLUDING REMARKS 

We have investigated the p-process nucleosynthesis in¬ 
side SSADs around a compact object of 1.4M0. Supernova 
explosion of 2O-4OM0 progenitors is probably responsible 
for the p-process nucleosynthesis inside SSADs. The chem¬ 
ical compositions of the accreting gas far from the central 
object have been assumed to be those of the oxygen/neon 
layers of the progenitor in the model of 10, 20, 30, and 
4 OM 0 stars. It is found that when to > 10® the temper¬ 
ature of the accreting gas attains to 2 — 3 x 10® K and 
appreciable nuclear reactions take place to produce the p- 
nuclei. We have also estimated the OPFs of the p-nuclei 
in the ejected material from the accretion disks via jets 
and/or winds to compare the abundances of the ejected 
matter with the solar abundances. The resultant p-nuclei 
profiles have similar features to those for Type H super¬ 
novae (Rayet at al. 1995). The light and intermediate 
mass p-nuclei such as ®^Mo, ^"‘Mo, ®®Ru, ®®Ru, and 

^®®La are underproduced. The amounts of p-nuclei and 
profiles of the OPFs ejected from SSADs depend on initial 
abundance distributions of s-process seeds. 

Moreover, for investigating the effects of proton and he¬ 
lium contamination via large-scale mixing in a supernova 
explosion of fallback material on the p-process nucleosyn¬ 
thesis, we have added some fractions of proton or helium 
to the initial composition for the cases of 10, 20, 30, and 
4 OM 0 progenitors and calculated the OPFs of the p-nuclei 
in the ejected material. It has been found that significant 
fractions of ^He inclusion cause small changes in the OPFs. 
Even for a small amount of proton included in fallback 
material, p-nuclei are synthesized through not only pho¬ 
todisintegrations but radiative proton captures. The most 
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prominent is an enhancement of ®^Mo. ^^®La is also over¬ 
produced for m = 10®. If fallback matter includes abun¬ 
dant proton, appreciable amounts of the p-nuclei, such as 
®^Mo, ®®Ru, ^^®Sn, and which are deficient in Type 

II supernovae, can be produced for m = 10®, though the 
overall profiles of the OPFs becomes worse compared with 
the solar ones. 

We have discussed the contribution from the p-process 
nucleosynthesis in SSADs to galactic chemical evolution 
of p-nuclei. We conclude that p-process in SSADs com¬ 
parably contributes to chemical evolution of p-nuclei as 
in supernovae if several percents of fallback matter are 
ejected from SSADs. In particular the p-process in SSADs 
may be important for the chemical evolution of ®^Mo and 
possibly compensates the underproduction in Type II su¬ 
pernovae. Contribution from SSADs to the chemical evo¬ 
lution of ®®Ru and ^®®La is unlikely to be significant. 

Numerical simulations of supernova-driven jets present 


the peculiar nucleosynthesis of the alpha-rich freeze out 
(Nagataki et al. 1997): The p-process will be affected by 
jets as SN1987A. Some gamma-ray bursts are likely to 
be caused by aspherical jet-induced supernova explosions 
(MacFadyen & Woosley 1999; MacFadyen et al. 2001). 
Numerical calculations indicate that in such explosions a 
SSAD is formed around a new-born compact object and a 
substantial amount of fallback (~ 0.2Mq) followed by mix¬ 
ing takes place (Fryer & Heger 2000; Hoflich, Khokhlov 
& Wang 2001). Hence some gamma-ray bursts possibly 
accompany with the synthesis of p-nuclei; Our SSAD p- 
process model would be responsible for the synthesis of 
p-nuclei in these bursts. 

We are grateful to Dr. Takashi Yoshida for his helpful 
processing nuclear data and useful discussion. We also 
thank the anonymous referee for useful comments. 
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Table 1 


Elements included in the nuclear reaction network. 


elements 

A 

elements 

elements 

A 


H 

1- 3 

Cu 

56- 71 

Ln 

121-139 

He 

3- 6 

Zn 

57- 74 

Ce 

126-142 

Li 

6- 8 

Ga 

60- 77 

Pr 

127-143 

Be 

7- 10 

Ge 

66- 76 

Nd 

132-150 

B 

8- 12 

As 

70- 77 

Pm 

133-151 

C 

11- 14 

Se 

71- 82 

Sm 

136-154 

N 

12- 15 

Br 

74- 83 

Eu 

137-155 

0 

14- 20 

Kr 

75- 86 

Gd 

140-160 

F 

17- 22 

Rb 

78- 87 

Tb 

143-161 

Ne 

17- 24 

Sr 

80- 88 

Dy 

146-164 

Na 

20- 27 

Y 

82- 92 

Ho 

149-165 

Mg 

20- 29 

Zr 

84- 96 

Er 

151-170 

Al 

22- 31 

Nb 

86- 97 

Tm 

152-171 

Si 

24- 34 

Mo 

88-100 

Yb 

154-176 

P 

27- 38 

Tc 

90-101 

Lu 

156-177 

S 

28- 42 

Ru 

92-104 

Hf 

158-180 

Cl 

31- 45 

Rh 

94-105 

Ta 

164-181 

Ar 

32- 48 

Pd 

96-110 

W 

166-186 

K 

35- 49 

Ag 

98-111 

Re 

168-187 

Ca 

36- 50 

Cd 

100-116 

Os 

170-192 

Sc 

39- 51 

In 

103-117 

Ir 

172-193 

Ti 

40- 53 

Sn 

105-124 

Pt 

174-198 

V 

43- 55 

Sb 

106-125 

Au 

178-200 

Cr 

44- 58 

Te 

108-130 

Hg 

180-204 

Mn 

46- 64 

I 

114-131 

Tl 

186-205 

Fe 

47- 65 

Xe 

116-136 

Pb 

190-208 

Co 

50- 66 

Cs 

117-137 

Bi 

192-209 

Ni 

51- 68 

Ba 

120-138 




Table 2 

Overproduction factor Fq ejected from SSADs for the 
PROGENITORS OF 10, 20, 30, AND 40 Mq AND FOR VARIOUS 
VALUES OF PROTON CONTAMINATION Xp. 


Mms 

(Mq) 




^p = 0 

Ap =0.01 

Ap = 0.1 

Ap = 0.3 

10 

18.15 

26.73 

21.40 

18.03 

20 

49.64 

161.1 

88.46 

89.03 

30 

102.0 

461.8 

259.1 

235.9 

40 

138.3 

698.5 

422.3 

332.2 


All Fq are calculated for the case of m = 10®. 
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Table 3 

Various masses of progenitors and fallback matter 


Mms’^ 


Afdest'^ 

WW95 model 

Mf‘1 Mrem" M^eed 

Mf‘1 

F99 model 

Mrem" Mseed^ 

15 

1.37 

1.55 

0.11 

1.43 

0 

0.2 

1.52 

0 

20 

1.67 

1.8 

0.32 

2.06 

0.26 




25 

1.48 

1.8 

0.29 

2.07 

0.27 

3.8 

5.58 

2.0 

30 



2.41 

4.24 





40 

2.36 

3.5 

8.36 

10.34 

4.84 

10.3 

11.01 

6.8 


Note. — All masses are in units of Mq. 

'^Progenitor mass on the main-sequence. 

’^Location of the base of oxygen burning shell from Hashimoto (1995). 

‘'Location where the peak temperature of supernova shock is equal to 2 xlO® K 
from Hashimoto (1995). 

‘^Fallback mass. 

®Remnant mass after fallback. 

^Estimated fallback mass with s-process seeds for p-nuclei. 



Fig. 1. — The initial abundances of models MIO (dashed line) M20 (solid line), and M40 (thick-solid line). The filled circles are the solar 
abundance (Anders Sz Grevesse 1989). 
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Fig. 2.— The abundance profiles of representative p-nuclei inside a SSAD for M20. The abscissa is the radius of the fallback disk in units 
of the Schwarzschild radius. The solid, thick-solid, dashed, thick-dashed, dotted, and thick-dotted lines represent ^“^Se, ^^Mo, ^^^La, ^^^Gd, 
^®^Ta, and respectively. 
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Fig. 3.— OPFs of 35 p-nuclei for M20 with m = 10®. Fq is equal to 49.64. 
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Fig. 4.— Same as figure 3, but for MIO. Fq is equal to 18.51. 
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Fig. 5.— Same as figure 3, but for M20P001. Fq is equal to 157.7. 
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Fig. 6.— OPFs of 35 p-nuclei for M20P01 with m = 10®. Fq is equal to 41.84 






Over production factor 



Fig. 7.— Same as figi 







